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Q\ ' Abstract 

The Faddeev-Yakubovsky equations have been solved in configuration space for the 
four nucleon system. The n+t elastic cross section has been calculated using several 
realistic interactions and a phenomenological three-nucleon force. Special attention 
is devoted to the description of the observed resonant structure at Ti a \, ~ 4 MeV. 
In the present state of the calculations, the realistic NN interactions are unable 
to describe the energy region considered, i.e. n+t S- and P- waves. Whereas the 
inclusion of a three-nucleon force improves the zero energy region, the disagreement 
in the resonant peak remains. This failure contrasts with the results of a pure 
phenomenological model. The origin of this disagreement as well as the possible 
t£) , issues are discussed. 
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1 Introduction 

*o 

The four nucleon continuum spectrum constitutes the present challenge of the 
few-body nuclear problem. Its interest lies not only in the natural progression 
that it represents towards the description of an increasing complexity but also 
in the richness of the A=4 nuclear chart itself. The problem implies, we believe, 
a qualitative jump with respect to the A=3 case. Indeed the failure of the usual 
NN interactions in describing the 3N system and the 4N bound state can be to 
a large extent overcome by the inclusion of three nucleon interactions (TNI) 
[1-5]. The success of this approach is mainly due to the strong correlation 
existing between the low energy scattering observables and the three nucleon 
binding energy [1,2,6], which in its turn is strongly correlated to the four 
nucleon one [7-9]. This means that the zero and low energy three and four 
nucleon problem depend on a few number of parameters which can be found 
in any existing TNI. By adjusting these parameters in order to reproduce the 
triton binding energy a global satisfactory description of the 3N system and 
the 4N bound state is reached [4,10,11]. The strength of the TNI depends on 
the underlying NN forces and their main effect is to rescale the 3N and 4N 
bound states energies. Their influence on the 3N continuum manifests itself 
only at relatively high energies [12]. 
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However, unlike the 3N case, the 4N continuum spectrum exhibits a rich va- 
riety of resonances and thresholds sufficiently far from the zero energy region 
that cannot be determined by the low energy properties. The simple compar- 
ison between the smooth behaviour observed in the n+d elastic cross section 
and the non trivial structure manifested in the n+t case illustrates this point 
well. It is therefore not clear that the approach followed until now, based on 
a good description of the deuteron and triton binding energies, could be suc- 
cessfully extended. A resonant state spreads over an energy region given by 
its width and fully displays there the internal dynamics of the 4N system. 
Its theoretical description provides a severe test in our understanding of the 
nuclear forces. 

Among the different scattering channels the n+t elastic case is the simplest 
one and so is the obliged first step before a more ambitious program can 
be undertaken. On the one hand, it is is a pure T=l channel free from the 
additional difficulties arising when dealing with Coulomb forces. On the other 
hand high quality measurement of the total n+t cross sections are available 
[13]. They show a resonant structure interpreted by the R-matrix analysis of 
Hale and collaborators as being generated by a family of resonances [14]. The 
aim of this work is to investigate whether the usual realistic NN interactions 
are able to account for these experimental facts. 

The results we will present in what follows have been obtained by solving 
the Faddeev-Yakubovsky (FY) equations in configuration space. The formal 
apparatus and details on the numerical methods can be found in [15,16]. We 
will restrict ourselves in this letter to a short presentation of the formalism in 
the case of four identical particles (next section) followed by the results and a 
conclusion. 



2 Formalism 

In the case of four identical fermions interacting via a pair-wise potential V, 
the FY equations result in a set of two integrodifferential equations, coupling 
two FY amplitudes K and H: 

(E-H - V)K = V [(P 23 + P 13 ) (-1 + P 34 ) K - (P 23 + P13) H] 
(E-H - V)H = V [(P 13 P 24 + P14P23) K + P 13 P 24 H] (1) 

, in which are the permutation operators. The wavefunction is given by: 

# = ^i +3 + ^ 2+2 
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^ 1+3 = [1 + e(Pi 3 + P23)] [1 + £(^14 + A 4 + P34)] K 

^2+2 = [1 + £(Pl3 + ^23 + ^14 + ^24) + ^13^24] # 



(2) 
(3) 



Each FY amplitude F=K,H is considered as a function of its own set of Jacobi 
coordinates x, y, z defined respectively by 



x K 

m 

z K 



r 2 -n x H = r 2 - n 

^ _ n±|±m) z H = V2 - ^) 



and expanded in angular variables for each coordinate according to 

(xyz\F)^J2 ^^-Y a (x : y : z). 
a xyz 



(4) 



The quantities F a are called regularized FY components and Y a are tripo- 
lar harmonics containing spin, isospin and angular momentum variables. The 
label a holds for the set of quantum numbers defined in a coupling scheme 
and includes the type of amplitude K or H. We have chosen the j-j coupling 
represented by 



K= {[(tih) Tx t 3 ] T3 U} T <g> {[(lx(siS 2 )* a ) jx (lyS3) jy ] j3 (IzSajj 



H: 



(f 1*2)7* (M^rjy ® I (lx(siS2)a x ) jx {ly{s^S A ) ay ) 



Jxy 



where s$ and ti are the spin and the isospin of the individual particles and 
T, J the isospin and total angular momentum of the four-body state. Each 
of the N c = Nk + Njj components in the expansion (4) is labeled by 12 
quantum numbers, restricted by the antisymmetry properties to the conditions 

(_)<T*+T«+J« = £f()rK (_y x +T x +l x = (_y y +T y +ly = £ for R 

The boundary conditions for a 1+3 scattering problem are implemented by 
imposing at large enough value of z the Dirichlet-type condition 



K(x,y,z) = t(x,y) 
H(x,y,z) = 

t(x, y) being the triton Faddeev amplitudes with quantum numbers (l x (siS2)a x )j x (h s ^)j y 
They ensure a solution which, e.g. for a relative n+t S-wave, behaves asymp- 
totically like 



J 3 
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K(x, y, z) ~ t(x, y) sin (qz + 5) 

where 5 is the n+t phaseshift and q, the conjugate momentum of the z-Jacobi 
coordinate in K-amplitudes, is related to the center of mass n+t kinetic energy 
T cm and the physical momentum k by T cm = |T lab = ^q 2 = \^k 2 



3 Results 

The results we will present concern the 4 He bound state and the n+t elastic 
scattering in the relative S- and P- waves. They have been obtained with the 
AV14 [17], NIJM-II and REID-93 [18] NN interactions. They turned out to 
provide very similar predictions in the observables we have considered and we 
have chosen to give in more detail only the results of AV14. They will also be 
compared to the predictions of the S-wave phenomenological MT I-III model 
published in [16]. 

For 4 He without Coulomb forces we obtain the binding energies (B) and r.m.s 
radius (f) displayed in Table 1. The NN interaction is limited to 1 S , 3 Si, Z D X 
partial waves and the number of FY components in expansion (4) is limited to 
N c = 5+5 and N c = 15+9, the same as those listed in [9]. The results provided 
by the different potentials are relatively close to each other but still far from 
the Coulomb corrected experimental value B = 29.0 MeV. For AV14 they are 
in good agreement with the momentum space calculations, e.g. B 5+5 = 23.36 
MeV and B 15+g = 23.77 MeV from [9], while more complete calculations give 
as converged values 24.2 MeV [33,20] and 24.7 MeV [32]. For the two remaining 
potentials we present here the first calculations. 
Table 1 

4 He binding energy B (MeV) and r.m.s. radius f (fm). 





N c = 


5 + 5 


N c = 


15 + 9 




B 


f 


B 


f 


AV14 


23.34 


1.56 


23.81 


1.54 


NIJM II 


23.39 


1.54 


23.86 


1.53 


REID 93 


23.65 


1.53 


24.12 


1.52 



In the scattering calculations we will expand the different J w n+t states in 
terms of the asymptotic hamiltonian eigenstates \L, S; J n > in which L is the 
n+t relative angular momentum and S = s n + s t the total spin. The following 
J 7 " states will be considered: 

|0+)= |0,0;0+) 
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= Ci+ 


|o, 


1 


l + ) + d 1+ 


2,1; 


1 + ) 


|0"> 




|1, 





0') 






in 


= Ci- 


|1, 


1 




1,0; 


I") 


|2"> 


= c 2 - 


|1, 


1 


2~) + d 2 - 


3,1; 


2-) 



(5) 



where the coefficients c and d, fixed by the dynamics, make them eigenstates 
of the S-matrix. These states are the only relevant ones for the low scattering 
energies considered. Moreover the coupling L <-> L + 2 between the different 
asymptotics in the J n = 1 + , 2~ states turned to be very small and has been 
neglected in the calculations presented below, which are thus restricted to 
the lowest angular momentum component in the asymptotic channel. For the 
remaining J n = 1~ coupled state we present the two corresponding eigenphase- 
shifts. 

3. 1 S-waves 

The phase-shifts for the first positive parity states (J 71 " = + , 1 + ) with the 
AV14 interaction are given in the first two columns of Table 2. They have 
been obtained with the NN potential acting only upon the 1 5'o, 3 Si, 3 Di partial 
waves and limiting the expansion (4) to l y , l z < 2 except in the K amplitudes 
where l z < 1. That makes for J w = + (l + ) a number of FY components equal 
to iV c = 19(48). This criterion corresponds to the N c = 15 + 9 calculations for 
the 4 He (see Table 1) and N c = 5 Faddeev amplitudes for the triton asymptotic 
state. This triton has a binding energy B$ = 7.43 MeV, to be compared with 
the fully converged value B 3 = 7.67 MeV obtained with N c = 34 [19]. 

Table 2 

AV14 n+t S- and P-wave phase-shifts (degrees) as a function of center of mass 
kinetic energy T cm (MeV). 



T 

± cm 


0+ 


1+ 


0+ (TNI) 


1+ (TNI) 


0" 


1- 


1- 


2~ 


0.05 


169.5 


170.8 


170.3 


171.4 


0.09 


0.15 


0.19 


0.11 


0.5 


147.9 


151.7 


150.1 


153.4 


2.2 


4.3 


5.9 


4.0 


1.0 


135.8 


140.8 


138.7 


143.1 


5.9 


10.5 


16.1 


10.9 


2.0 


120.5 


126.9 


124.1 


129.7 


13.8 


20.6 


35.9 


25.0 


3.0 


110.2 


117.2 


114.0 


120.3 


21.2 


26.6 


47.8 


35.4 


4.0 


102.3 


109.8 


106.3 


113.0 


27.4 


30.1 


53.5 


41.9 


5.0 


95.5 


103.7 


99.9 


106.9 


32.3 


32.1 


56.3 


45.8 


6.0 






94.3 


101.6 


34.8 


35.1 


59.6 


48.7 



The low energy parameters are given in Table 3 for several potentials. The 
scattering lengths have been calculated directly at zero energy whereas the 
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effective range is extracted from fitting the results of Table 2. One can see, as 
in the 4 He bound state, very similar results for the potentials considered. Some 
dependence on the maximum value of l z is observed. The values for AV14 are 
in agreement with those obtained with a variational approach a(0 + ) = 4.32 
fm and a(l+) = 3.80 fm [20]. 

Table 3 

Singlet, triplet and coherent n+t scattering lengths a(0 + ), a(l + ), a c and effective 
range tq (fm). The zero energy cross section ctq is in fm 2 . 





lz 


= 






h = 


0,1 








a(0+) 


a(l+) 


a(0+) 


r (0 + ) 


a(l+) 


ro(l + ) 


a c 




NIJM II 


4.34 




4.31 




3.76 




3.90 


191.6 


AV14 


4.34 


3.87 


4.31 


2.08 


3.79 


1.76 


3.92 


193.7 


AV14+TNI 






3.99 


1.95 


3.53 


1.71 


3.65 


167.5 



The S-wave n+t elastic cross section is displayed in Figure 1 for AV14 (dotted 
thick line) and compared to the experimental values taken from [13]. We have 
included for comparison the results obtained with MT I-III potential (solid 
thin line) [16]. The separated singlet (long-dashed) and triplet (short-dashed) 
contributions are also plotted. The disagreement with the experimental data 
is dramatic although expected due to the strong correlation between the 4N 
binding energy (far from its experimental value) and the N+3N scattering 
length. This disagreement can be overcome as in the 3N case by including a 
three nucleon interaction. Several TNI models exist in the literature [21-25]. 
In all of them, some parameters have had to be adjusted in order to reproduce 
the experimental 3N, and eventually 4N, binding energy. In view of that we 
have used the phenomenological model having the form 

W(p) = W r W a , p = J x 2 + y 2 (6) 



The parameters have been chosen to reproduce the experimental triton binding 
energy B 3 = 8.48 MeV. With the parameter set W r = 500 MeV, W a = 174 
MeV, p: = 2.0 fm, one has also a satisfactory agreement for the 4N binding 
energy = 29.0 MeV and the following n+t low energy parameters: ao = 4.0 
fm, ai = 3.53 fm and er(0) = 167.5 fm 2 . We notice that the parameter set 
used is far from being unique and that the four-body calculations have been 
performed with N c = 5 + 5 FY components. The corresponding n+t S-wave 
cross section has also been calculated and is plotted in Figure 1 (solid thick 
line). The zero energy cross section is now in agreement with the experimental 
points. We notice also that the values of the scattering lengths thus obtained 
are close to those provided by the variational CHH method, ao = 4.12 fm, 
a\ = 3.59 fm and cr(0) = 172.4 fm 2 , adding to the AV14 NN potential the more 
sophisticated Urbana TNI-models [20]. This seems to prove that, at least for 
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the observables under investigation, the essential role of a three nucleon force 
is just to provide the right three {B 3 ) and four body (B4) binding energies 
and consequently ensure acceptable n+t low energy parameters. 




Fig. 1. n+t S-wave cross section with MT I-III, AV14, and AV14+TNI interactions. 

The experimental situation concerning the n+t scattering lengths is not very 
well established (see Table 4). Only the zero energy cross section a exp {Q) = 
1.70 ± 0.03 b [13] seems to be a firm result. The coherent scattering length 
a c from [26] is based on the assumption ai/a = 0.92 ± 0.04 and the deduced 
values for a® and a\ are thus not purely experimental results. The more recent 
measurement [27] gives a c = 3.59 ± 0.02, a value just compatible with the 
preceding one, but the two solution sets inferred for the scattering lengths 
are too far from any theoretical prediction to take them without some care. 
The values given in [28] are in fact extracted from p- 3 He data via a Coulomb 
corrected R-matrix analysis. The results of our calculations together with the 
MT I-III predictions tends to favour values close to a® = 4.0 — 4.1 fm, a\ = 
3.53 — 3.60 fm and a c = 3.65 — 3.70, i.e. the older experimental result of [26]. It 
would be however of a great interest to have a definite experimental conclusion 
on the a c value. 
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Table 4 

Experimental coherent (a c ), singlet (cto) and triplet (ai) n+t scattering lengths (fm). 



a c = (3ai + a )/4 


a 


ai 


ref. 


3.68+0.05 


3.91±0.12 


3.60±0.1 


[26] 


3.59+0.02 


4.98+0.29 


3.13±0.11 


[27]-I 




2.10±0.31 


4.05+0.09 


[27]-II 


3.607+0.017 


4.45+0.10 


3.32+0.02 


[28] 



3.2 P-waves 



The phase-shifts for the J 71 " = 0~, 2~ states, are given in Table 2. The two 
J w = 1~ columns are the eigenphase-shifts corresponding to the two coupled 
asymptotics states in (5). 

3.0 



2.5 



2.0 



1.5 



1.0 



0.5 



0.0 




- P-waves (AV14) 

- S+P (AV14) 

- S+P (AV14+TNI) 
+ Experiment 



/ 



0.1 



1.0 



T lab (MeV) 



Fig. 2. n+t cross section with AV14 with and without three nucleon interaction. 

The resulting n+t elastic cross section are displayed in Figure 2. It appears 
that the results obtained with AV14 (dashed line) are unable to account for 
the experimental data. By removing the S-wave part in the cross section one 
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finds that ~ 1/3 of the P-wave contribution is missed in the calculations. This 
is in contrast with the MT I-III results, which provides an almost perfect de- 
scription of the resonant region [16]. It is worth noticing that in the presence of 
the phenomenological three-nucleon forces (6) this quantitative disagreement 
remains and is even worse (solid line). Whereas the TNI sensibly rescales the 
zero-energy cross section, its doesn't have any significant effect on the P-wave 
phase-shifts. The differences between both cross sections, visible in Figure 2, 
are essentially due to the n+t S-wave contribution. The small influence of the 
three nucleon interactions in the low energy scattering observables seems to be 
a general property of the existing models [12] rather than a consequence of the 
simplified TNI we have considered in (6). A similar disagreement can be ob- 
served in the differential cross section at neutron laboratory energy Ti^ = 3.5 
MeV taken from [29] and displayed in Figure 3. 



500 



400 



-I 300 



a 

-a 

200 



100 







Experiment 
AV14 
AV14+TNI 








30 



60 



90 

9 (deg) 



120 



150 



180 



Fig. 3. The n+t differential cross section at Ti a b = 3.5 MeV obtained with AV14 
alone and including three-nucleon interaction are compared to experimental data. 

The comparison between the different P-waves contributions given by MT I-III 
and AV14 potentials is done in Figure 4 and shows that: 

(i) the degeneracy between the J 71 = 0~ , 1 _ , 2~ states existing in the MT I-III 
potential no longer holds in the case of AV14 as a consequence of its tensor 

— # 

and / • s terms. This splitting broadens the sum of these three contributions, 
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which present their maximum separately and results in a lower value for the 
cross section. 

(ii) there is additionally a real lack in their separated contributions, specially 
in the 0~ and the statistically enhanced 2~ partial wave. 

1.5 



1.0 



0.5 



0.0 




T Iab (MeV) 



Fig. 4. Comparison of the cross sections obtained with MT I-III and AV14 for each 
partial wave. 



A possible influence of P-wave NN interaction has been checked in [30]. The 
1 Pi, 3 Po, 3 Pi NN partial waves were added to 1 S , 3 Si 3 D\, thus completing 
the shell, but the corrections were found to go in the wrong direction. 
Since, we have included the 3 P2 — 3 F2 components and we present in Table 5 
the phase-shifts at T cm = 3 MeV. The corresponding elastic cross sections, 
togheter with the MT I-III and NIJM II predictions, are given in Table 6 and 
plotted in Figure 5. Even if the relative weights of the different partial waves 
are modified, the disagreement noticed so far still remains. The inclusion of the 
different NN P-waves interactions has no global effect in the cross section. This 
conclusion is in agreement with the extensive study done for 3N continuum 
case (see e.g. Fig. 6 in [10]) in which the was shown to give satisfactory 
results at the considered low energies. 
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Another reason for such a disagreement, suggested by the success of the MT 
I-III model, could be found in the NN S-wave tensor force itself. In partic- 
ular we managed to increase the n+t P-waves phase-shifts and consequently 
the total cross section by weakening the 3 Si — 3 Di coupling in benefit of the 
3 S\ potential, keeping the same triton binding energy all along. The tensor 
force, firmly established in the deuteron and the NN observables, seems over- 
estimated in the more compact systems and responsible for the underbinding 
problem in the 3N and 4N bound states. A possible solution could be found 
by suppressing its the short range part. 

Table 5 

n+t P-waves phase-shifts (degrees) at T cm =3 MeV, with different truncations of the 
NN interaction. 



NN partial waves 


o- 


1- 


1- 


2" 


1 5o, 3 Si, 3 D\ 


21.2 


26.6 


47.8 


35.4 


3 5 l5 3 D, + ip l5 3 P , 3 Pi 


24.1 


20.0 


37.5 


32.7 


3 5 l5 3 D 1 + 1 P 1 , 3 P , 3 Pi, 3 P 2 , 3 F 2 


25.1 


22.8 


38.0 


41.4 



Table 6 

n+t elastic cross section a (barns) at T cm =3 MeV for different potentials. 

a(b) 



exp. 2.40 

MT I-III 2.40 

AV14 (no TNI, no P-waves) 2.12 

AV14 (TNI, no P-waves) 2.06 

AV14 (TNI, j x < 2) 1.78 

AV14 (TNI, all P-waves) 2.03 



NIJM II (no TNI, no P-waves) 2.09 



Under the restrictions in which they have been obtained, i.e. limited number 
of amplitudes in expansion (4), the preceding results indicate that the AV14 
interaction, limited to NN S- and P-waves, cannot reproduce the n+t elastic 
cross section. The calculations performed with other realistic local potentials 
lead to very similar results for the scattering lengths (Table 3) as well as for 
the phase-shifts (Table 7). The qualitative behaviour of AV14 seems thus to be 
a general feature of the modern realistic NN potentials at least in their local 
form. On the other hand the recent results obtained by Fonseca [31] reach 
similar conclusions in a totally independent way and extend them extensible 
to the non local Bonn B potential. 
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O AVI 4 (TNI, no P- waves) 

AAV14(TNI,j x <2) 

V AVI 4 (TNI, all P- waves) 

A NIJM II (no TNI, no P-waves) 



1.0 



T,. t (MeV) 



+ 
+ 



lab 



Fig. 5. Summary of different calculations for the n+t elastic cross section at 4 MeV 
laboratory incident energy. 

Table 7 

Comparison of the n+t P-waves phase-shifts (degrees) at T cm =3 MeV with AV14 
and NIJM II. The NN interaction is restricted to its 1 Sq, 3 S\, 3 Z?i components. 



NN potential 


O - 


l - 


1" 


2 - 


AV14 


21.2 


26.6 


47.8 


35.4 


NIJM II 


20.6 


26.0 


46.0 


34.1 



4 Conclusion 



The solutions of Faddeev-Yakubovsky equations in configuration space for 
the low energy neutron-triton scattering have been obtained by using several 
realistic NN interactions. 

The results obtained for the n+t scattering length and elastic cross section 
are compared to the existing data. The interactions considered, at least when 
limited to the NN S- and P-waves, fail in describing both the zero energy cross 
section and the n+t resonant structure observed at neutron Ti ab ~ 3.5 MeV. 
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A phenomenological three-nucleon interaction has been included. It provides 
overall agreement for 3 H and 4 He binding energies as well as for the n+t S-wave 
scattering region. Its effect is however negligible in the resonant region and the 
disagreement remains. Including NN P-wave interactions doesn't substantially 
modify this result. 

On another hand the success obtained by using simple S-wave spin dependent 
potentials [16] leads us to think that the combined effect of Vnn + Vnnn + ■ ■ ■ 
forces results in a simplified effective interaction, as happens in the nuclear 
structure calculations, and shows how little effectives the existing OBEP mod- 
els are in describing the few nucleon systems. 

These results need to be confirmed and the convergence in expansion (4) 
pushed further. In their present state they suggest a basic inadequacy of the 
existing NN interactions in describing the n+t resonant states, a problem 
which can hardly be corrected in terms of three nucleon forces. 
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